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SUMMARY 

The important role of Mathematics in the scientific development of human- 
ity is determined by the possibility of elaborating mathematical models of the 
objects studied in the different spheres of science. In other terms, it means, 
describing by means of the rigorous language of mathematics the properties 
and relationships of these objects, that in turn it allows the application of the 
rigorous and precise techniques of this science to the problems related with the 
study of the real objects facilitating this way the adquisition of conclusions on 
these and even to control their behavior. 

The application of these models and methods of mathematics have been 
increased significantly from the second half of the XX century, due to: the 
appearance and development of the calculation techniques, the introduction of 
the systemic focus in the sciences and the development of new branches that it 
constitutes the Modern Applied Mathematics. 

The fundamental objective of mathematics in the formation of specialists 
in the different branches of the science and technical areas is to establish that 
they are capable of applying the mathematical modelation and the calculation 
experiment, like new technology of the scientific work in all the spheres of the 
professional activity. To achieve this objective it is necessary to develop in the 
students abilities for: the application of logical forms of thinking, the application 
of the mathematical concepts in the description of the objects and relationships 
among the real objects and for the elaboration, selection and application of 
the most efficient algorithms for the computational solution of the outlined 
problems. [1] 

The identification theory and realization of the dynamic systems is a medullary 
aspect in the modern control theory that consists fundamentally in that, start- 
ing from the knowledge of the behavior entrance-exit, obtained experimentally 
in the case of the identification, or given previously in the case of the realization, 
to build a state model that carries out this behavior. This content is not gen- 
erally treated in the pre-graduate courses, for the systems of multiple entrances 
and multiple exits. 

In this work it is demonstrated that the identification theory and realization 
of the lineal dynamic systems can be imparted in the technical careers starting 
from the results of Mathematical Analysis and Lineal Algebra received by the 
students in pre-graduate studies, without necessity of adding new contents in 
the programs of this subjects 

We propose a new form of imparting the theory of identification and real- 
ization of the lineal dynamic systems based on the intuition and the physical 
interpretation of the concepts. 

INTRODUCTION 
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The daily work of mathematicians necessarily consists in obtaining new the- 
orems, the discovery of new connections among concepts already well-known 
and the introduction of new concepts. However, this is not the only objective 
of mathematical investigation, the academic A. N. Kolmogorov also considered 
as important objectives of mathematical investigation the following : 

1-Carry the modern logical foundations of mathematics to such a level that 
they can be exposed the scholars from 14 to 15 years. 

2 - To destroy the divergences among the rigorous methods of the pure 
mathematicians and the non rigorous focuses of the mathematical reasonings 
used by the applied mathematicians, the physicists and the technologists. [2] 

With relationship to these objectives the mathematicians have the task of 
elaborating the logical foundations of the irreproachable mathematical intu- 
itions, in which the rigorous calculation methods are based, arosen from the 
alive intuition of the physicists and the technologists that have been validated 
in the practice. But these foundations lead in occasions to constructions so 
sophisticated and complex that make the mathematicians proud, but that the 
physicists and engineers would not have the possibility to assimilate if the nec- 
essary effort was not carried out to make it accessible to them 

This work meets the first objective of the investigation outlined by Kolmol- 
gorov and has as a result: 

To develop the identification theory and realization of the lineal dynamic 
systems starting from the results of Mathematical Analysis and Lineal Algebra 
received by the students in pre-graduate studies with a methodological focus 
that facilitates their assimilation for the students starting from the pre-grad 
level. 

DEVELOPMENT 

In this work a methodology is presented for the exposing of the problem of 
obtaining a state model of a dynamic system, starting from the knowledge of its 
entry- exit bahaviour, given by a transfer matrix, so that this is understandable 
by the students from the pre-grad level whose mathematical foundation is based 
on results of Functional Analysis and of the theory of matrixes whose elements 
are rational functions. With this methodology it is avoided to use sophisticated 
and complex mathematical theories. 

The structural unit and logics of the exhibition is guaranteed by means of the 
application of analysis procedures, synthesis, comparison exemplification and 
of the intuition starting from the physical interpretation of the mathematical 
concepts and the mathematical description of the physical problems. 

MODELLING OF SISO SYSTEMS . 

Starting from the analysis of different cases already well-known, as the one- 
dimensional movement of a body under the action of an external force, the 
circulation of the current in a circuit R-L, due to the application of a voltage 
source in the net, studied inside the Physics, etc. will establish that the students 
can arrive at the conclusion that the modclation of many technical problems 
leads to the position of a system with a simple entrance and a simple exit 
(SISO) 

Where: 
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u(t): are the variables of entry of the system that contain the information 
on the environment of what is disposed and the impacts in the behavior of the 
study object. 

y(t): are the variables of exit of the system that contain the information 
which is able to be captured through the sensors. 

x(t): variables of state of the system that contain the minimum quantity of 
information on the system that is enough to predict their future behavior. 

If the physical laws that govern the behavior of the real system are known 
and also this is does not vary with regard to the adjournments in the time and 
it satisfies the overlapping principle, that is to say it is stationary and lineal, 
and it is also initially in rest, the mathematical pattern is obtained: 

(D n + a 1 D n - 1 + ... + a n _ 1 D + a n ) X (t) = u(t) , X(0) = (1) 

y(t) = (b Q + b 1 D + ... + b k D k ) X(t) , k < n (2) 

Through the application of Laplace Transformation , the results obtained 

are: 

Y(S) = G(S) U(S) (3) 
where: 

G(S) = „„ + 6 i 1 i + "- + bk s * ( 4 ) 

it is denominated function of transfer of the system (which can be interpreted 
as the transformation of Laplace of the answer to the unitary impulse) 

Using the transformation of coordinated , x J_1 (t) = jq (t) i=l,- ■ • , n, the 
equations (1) and (2) can be expressed by means of the normal system of order 
n 

X = AX + BU 
Y=CX (5) 
where: 
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being I„_ithe unitary womb of size (n -1), X G R n , A G M (n ,n), 

B € M(n,l) and C S M(l,n). The womb A it is known in the literature like 
main partner of the characteristic polynomial of the equation (1) 

We can observe that the transfer function G(S) it is a rational function where 
the degree of the polynomial of the denominator coincides with the number of 
variables that intervene in the system and the polynomial of the numerator is 
associated to the exit of the system .Comparing (4) and (5) it is easy to establish 
a correspondence biunivoca among the pattern entry/exit given by the transfer 
matrix and the state model given by the tern (A, B, C) and pass from one to 
another without using the equations (1) and (2). 

This type of model is denominated "state model" and the variable x^ ; i 
= 1,- • • ,n ; are denominated state variables. 

MODELLING OF MIMO SYSTEMS 

Many of the problems that are presented in practicals respond to a model 
of multiple entrances and multiple exits (MIMO). We will Consider that these 



3 



systems satisfy the overlapping principle and they are invariable with regard to 
the translation in time, that is to say that they are lineal and stationary. We 
will also suppose that it demands the maximum degree of pairing between the 
entrances and the exits. Then for the linealidad of the system it can planted 
that: 

y» = Un + lh2 + ■■■ + Vij + ■■■ + Vir (6) 

where y^ represents the effect of the entry (j-esima) in the exit (i-esima), 

this jeans that it is the signal of the exit (i-esima), if only the entry (j-esima) 

was not null 

Uj ^ 0; u\ = u 2 = ... = uj-i = ... = u n = 0) (7) 

The same as in the case of the systems SISO where initially an analysis of 
the system in rest was made, that is to say for those in u(t)=0, y(0)=0 for t<0 

If we apply the Transformation of Laplace to the entrances and the exits of 
(7) we can obtain the transfer function that relates the exit i-esima with the 
entrance j-esima 

G v = vm ^ YiiW = Ga (S) Uj(S) (8) 

After if we apply Transformation of Laplace to (6) and we use (8) we will 
obtain 

Yi(S) = Ed, (S) Uj(S); i = l,..,m 

3=1 

whose matrical form is: 
Y(S) = G(S) U(S) (9) 
where: 

G(S) = [G i j (S) ] ( i =l,...,m ; j =1,..., r ) is a matrice of size m x r,Y(S) 
G R m y U(S) e R n 

Notice that equation (9) is the multidimensional analogy of equation (3) 
therefore we can say (including the multidimensional and the unidimensional 
case) that the Laplace transformation of the exit vector in a lineal system is a 
lineal function of the Laplace transformation of the entrance vector, defined by 
a matrice whose elements are functions of the variable (s) and that it is denom- 
inated transfer matrice. The physical interpretation of this matrice can also 
be made intending to apply to each entrance for shift a very narrow and sharp 
pulse, of unitary area, that ideally is equal to apply Ufc(t)=5jy5(t) in the ex- 
periment j-esimo, k,j=l,...,r, in each experiment the exits are measured and the 
observations Jij{t), for i=l,,m and j=l,...,r. The matrice [y»j(t)] is the answer 
matrice to the unitary impulse and the Laplace transformation [L _1 {yjj(t)}] = 
G(s), is the transfer matrice. 

To continue , the following problem is planted: ^How to obtain the transfer 
matrice starting from the state model? 

Be a state MIMO model: 

X = AX + BU 

Y=CX 

Applying the Laplace Transformation, we obtain: 
S X(S) - A X(S) + B U(S) 
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Y(S) = C X(S) 

Making X (S) the subject of the formula in the first equation and substituting 
in the second equation we obtain: 
Y(S) = C(SI- A)- 1 BU(S) 
The transfer matrices is given by: 
G(S) = C(SI-A)- 1 B 

This means that the solution of the problem outlined only requires algebraic 
calculations of inversion and multiplication of matrices that can be solved with 
the help of MatLab. 

REALIZATION OF STATIONARY LINEAL SYSTEMS 

Let us analyze the inverse problem: Given the transfer matrice G(S) = Gi 
j(S) (i=l,...,m; j=l,...,n) obtained using the formula experimentally (9) and 
carrying out the measurment that were indicated when giving the geometric 
interpretation of the transfer matrice. 

,/How to obtain the state pattern? 

This inverse problem that one knows with the name of realization problem, 
has great importance from a practical point of view and its physical interpreta- 
tion would be the following : 

Given a real system, with consistent uncertainty in that we don't know the 
physical laws that govern their development in time neither even define their 
state variables, we only have the possibility to excite the system with certain 
entrance signs and through the measuring of certain exit signs inorder to analyze 
their behavior . 

We will use an inductive method, that is to say, through the analysis of 
particular cases we will try to infer a general methodology 

Be the system: has two entrances and an exit, this means that the inves- 
tigator can apply two types of entrance signs, for example (an electric field 
and a magnetic field, or pressure and temperature) and it can measure an exit. 
Therefore the transfer matrice has the type: 

[Gll(S) G12(S)] 

By means of the formula (9) and the experiments described in the physical 
interpretation of the transfer function, it is concluded: 

ler experiment: m=#(t), U2(t) =0, y(t) = Exp(-t), then Gn (t) = jj^j = 

l 

s+i 

2do experiment: Ui(t)= 0, u 2 =S(t), y(t)=Cosh(t), then G12 (i) = ^(S) = 

S 

Then 

Comparing the results obtained in (4) and (5) we arrive for analogy to the 
following submodels: 

The corresponding to Gn represented by the tern (An, Bn,Cn) where: 
An = (-1); Bu = (1) and Cn = (1) 

And the corresponding to G12 represented by the tern (A12, B12, C12) where: 
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Al2 = 



and C ( 1 1 ) 



\ ) ' Sl2 = ( 1 
Then the resulting system is: 

ii = — Xi + Ml 

£2 = 23 

X3 = X 2 + U2 

yi=xi+ x 3 

whose state model is defined by the tern (A, B, C) where: 



A=| 




and C = ( 1 1 ) 



/ 1 
B = 

\o 1 

Inorder to synthesize and to generalize the obtained methodology it will be 
requested the students make an elaboration of an algorithm like the following 
one: 

Entry of data: 

6 (M) + 6 (i,J) g + ... +b (.i : i)gHj 



To receive G = [G ij ] (i = l„m ; j = l„n) ; Gy - - a (.. 3 ) s ^-i + . +a (i ;J ) 

Construction of sub models: 
For each j = l,...,r 
For each i = l,...,m 
To construct A i j eM ( n i j, n i j ) 

Make the last line equal to: 



-a; 



Aid) 



( 



Make the first column equal to: 



Complete the matrice A i j with the unitary matrice l nij -\ 
Define 















Define 

Increase i 
Increase j 

Make A equal to the cellular diagonal matrice with cells Ay , ( i=l,...,m ; j 
l,...,r) 
Make 



6 



6 (2,D 
b (n,l) 









make 



B= 








6(1.2) 
6(2,2) 



C 



( C W 


V o 





c (2,l) 












6(V 

6(2,r 
6(m,r) y 





r (m,l) 



c (l,2) 










c (2,2) 








o(m,2) 



c (l,r) 










c (2,r) 



This methodology will be able to be applied always and provided that they 
are rational fractions, that which is presented in most of the practical examples. 
If some element of G(s) it is not a rational fraction, then the system can not be 
modeled with equations of the type (4). 

MODEL OF MINIMAL STATE 

Until this moment we have obtained a solution of the problem for the cases in 
which the elements of the transfer matrice are rational fractions characteristic 
of the variable s , but can we guarantee that the solution is unique? and in 
negative case ;Is that the most simple solution? We will try to respond these 
questions by means of examples, avoiding the demonstration of a mathematical 
theorem where it can be concluded that the solution of the equation matricial 

C ( SI- A y 1 B = G(S) is not unique 

Above all it is defined when two lineal dynamic systems are equivalent, which 
is refered to as the pair 

(t, x) e (R, X ) where X is the state space, a phase and R x X phase space 

Definition: Two lineal dynamic systems with state vectors x y x are alge- 
braically equivalent whenever their numeric phase vectors are related for all t 
like (t,x)= (t, T x) where T is a non singular matrice of size n. [3] 

In the cases that we are analyzing the system of differential equations is 
lineal and stationary. Be this the state model abtained 

X = AX + BU 

Y=CX (A) 

Making a base change in the state space by means of the non singular trans- 
formation T, Z=TX we obtain a system algebraically equivalent to (A) 
Z = AX + BU 
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Y=C Z (B) 
where 

A= T A T _1 , B = TB , C = CT' 1 

It is demonstrated that two systems algebraically equivalent have the same 
transfer matrice and it can be concluded ternas (A, B,C) y ( A, B, C) and are 
solutions different from G = C ( S I - A ) _1 and as a consequence this equation 
has infinite solutions, since starting from any main non singular matrice T we 
can obtain the corresponding ones ( A, B, C) . 

This means that by means of the algebraic equivalence we obtain another 
state model, and therefore this is not unique. For these solutions the state 
space has the same dimension. KALMAN proved that there are states with 
different dimensions for the same transfer matrice and it outlined the problem 
of extracting the pattern of state of minimum dimension, that is to say, given 
a transfer matrice G(s) to obtain the state pattern with smaller number of 
variables of possible state. [3] 

To obtain the minimal state it is necessary to define what we call completely 
controllable system and system totally observable. 

Definition : A lineal dynamic system is completely controllable if it is not 
algebraically equivalent to a system of the type: 

Z\ = An zi + A12 z 2 + Bi u 

Z2= A22 Z2 

Y = ci zi + c 2 z 2 

where zi y z 2 are vectors of ni y n 2 = n - ni components respectively 
In other words, if the lineal system (A) is algebraically equivalent to a system 
that has for state model terna ( A, B, C) such that: 

H A £) ) » e -*« 

The system (A) is not completely controllable, that is to say, it is possible 
to find a system of coordinates in which the variable Zj is separated in two 
groups zi = [z 1 , z 2 , • • • , z™ 1 ) , y z 2 (z ni+1 , • • • , z n ), which are denominated 
controllable and not controllable respectively. 

The application of the definition is annoying since if the initial system has 
the form of ( A, B, C) it is evident that it is controllable, but if it doesn't 
have, we can not affirm nothing since there are infinite possibilities of algebraic 
transformations. For such a reason we will apply without demonstration a 
characterization based on the concept of range of a matrice. We could have used 
the variant of defining the concept however by means of the characterization that 
has possible algorithms, it is a very abstract definition. The form that we use 
gives the possibility to give an intuitive idea, a physical interpretation of the 
concept and for this reason we have selected it. 



The caracterization is the following: 
In the system: 
X = AX + BU 

Y=CX (A) where AeM (n, n ) 
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We consider that B is a matrice of the size n x r and M c = [ B AB- • • A™ 1 
B ] a matrice of n rows and n x r coloums called controllable matrice. 

" The condition necesary and sufficient for the system (A) to be totally con- 
trollableis that 

R(M c ) = n " [4]. 

To obtain the R(M c) an algorithm is deduced to find the k ( k < n ) columns 
lineally independent of M c without working jointly with M c columns. 
To reciever A y B 

To define H, a group of at most n vectors R™, initially empty 
for a= 1 to a= r choose b a e B 
make j = 

while r [ b a Ab a A j b a ] = j+1 

Include the vector coloumn PJb a in H 

Increase j 

Increase a 

Return H 

End 

Definition : The space generated by all the columns lineally independent of 
the matrice M c is called controllable space and we will denote it as E c . 
Property: E c is an invariant [5] with relationship to the system 

X = AX + BU 
Y=CX 

/How to separate the completely controllable and not completely controllable 

spaces? 

If E c = < fi • ■ ■ f k > being { fi • • ■ f &} a base of E c and let us complete with 
the canonical vectors { ei,- • ■ e„ } a base of En. If the lineal transformation 
T=[fi ••• fk e k+i ■ ■ ■ e n] of E n such that 

Z = T X = xi fj H + x k f fe+ x k +i Cfe+i H + x„ e „ 

It is demonstrated that by means of the matrice T a system is obtained 
algebraically equivalent to (A) whose state model is the terna ( A, B, C) where: 

H A t) iS =( B o") ?/ C whatever. 
It means, 

Zi = An zi + A12 z 2 + Bi u 

Z2= A.22 Z2 

Y = Cz 

We will now analyze the dual concept with regard to the controllable concept. 
Definition : A lineal dynamic system is totally observable if this it is not 
algebraically equivalent to a system of the type: 
Z t = An zi + Bi u 
Z 2 = A 2 i zi+A 2 2 z 2 + B 2 u 

Y =c\Zi 
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where zi y Z2 are vectors of niand 112 = n- ni components respectively. [3] 
In other words, if the lineal system (A) is algebraically equivalent to a system 
that has for state model the tern where: 



A=( . 11 9 ] ; C — (ci 0) and B whatever 
\ A 21 A22 J 



The system (A) it is not totally observable, that is to say, it is possible 
to find a system of coordinates in which the variable Zj is separated in two 
groups zi = (z 1 , z 2 , ■ ■ ■ , z ni ) , and z 2 • • • , z n ) which are denominated 

observables and non observables respectively. 

The same as in the study of the control the application of the definition is 
annoying and for such a reason we apply the following characterization: 
In the system: 
X = AX + BU 

Y=CX (A) where A e M (n,n ) 

Consider that C is a matrice of the size m x n and 
( C \ 



CA 



A matrice of n m rows and n coloumns called matrice of observability 
"The condition necesary and suficient for the system (A) to be totally ob- 
servable is R (M ) = n " [4] 

As transposing a matrice its range don't vary, it have: 

r(M ) = r ( C" A*C* • • • {A^^C 1 ) 
Note: The system: 
X = AX + BU 
Y=CX 

It is totally observable if and alony if X= A* X+ C*U is completely control- 
lable, in this sense it is to say that both concepts are dual. Consequently the R 
(M ) is calculated equally to R( M c )■ 

In this case the space generated by the lineally independent rows of Mo is 
denominated observable space and we will denote it E o and their complement 
ortogonal is denominated non observable space (E no ) 

How to seperate the sub-spaces E no y E o ? 

Let us suppose that Ec = [gi ■ ■ ■ g n -ki] and E^o = [vi ■ ■ ■ v^] , if wc 
complete the base in the non observabile sub-space to a base in the whole space 
wc obtain a matrice Ti. 

If the lineal transformation Ti = [vi • • • Vk et+i • • • e„] of E n such that 

Z = Ti X = vixi +• • • + v fcl x fel + v fel+ i c fel+ i +• • • + v„ e „. 

Let us demonstrate that by means of the matrice Ti a system is obtained 
that is algebraically equivalent to (A) whose state model is the tern ( A, B, C). 
where: 
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B whatever and C — [0 C 2 ] , or if, 



j_( Mi A 12 

V ^22 

Zi = An zi + A12 z 2 + Bi u 
Z 2 = A 2 2 z 2 + B 2 u 
Y = C 2 z 2 
Dcmostration: 

In effect: E no = { z : z 2 = } where Z = 

1" $")> B -(f 

hi A22 I \ r>2 



Zi 

Z2 



If A= ( . . J ; B = I ) and C = [ci c 2 ] , therefore the system 



is: 

Z\ = An zi + A12 z 2 
Z2 = A 21 z i+A 22 z 2 
Y = ci zi + c 2 z 2 

We choose z £ E no where Z = ^ Z q° ^therefore z(t) GEc; Z(i) = ^ q 1 

belong to E no then, ci zio =0 

Making z equal to the vectors canonicals of E no in z} ci+z 2 C2H l- z io c fei = 



we obtain: 

ci = c 2 = • • • = c kl = 
and C = [0 C 2 ] 

As z 2 = =>■ z 2 = => A 2 i zio = => A 2 i = 0, as z 10 can be whatever 
arbituary in the sub-space of non-controllability. Then 

A=[ J \ 12 I obtained the desired 



A 22 

Using analysis procedures and synthesis we can outline the following method- 
ology: 

ALGORITHM FOR THE CANONICAL DECOMPOSITION OF A SISTEMA MIMO, LINEAL AND SI 

Inicial data: A G M (n, n), B e M (n, r), C G M (m, n). 

Start 

1- Calculate: rank [ B AB • • • A" -1 B ] = k and extract a base from the 
sub-space E c generated by the coloumns of this matrice: 

E c = (ei • • • e k ) , dim E c = k -< n, E c C n n ,E G M (n, k), 
E = [ ei • • • e k ] 

2- Calculate rank( C" A*C* • • • (A t ) n ~ 1 C t ) = p and extract a base 
from the sub-space E o generated by the columns of this matrice: 

E o =(/! ■■■/„), dim E o = p < n, E„ C E F G M (n, p), F = [ fi ■■■ f p ] 

3- If E no = E|. calculate a base from the sub-space E c n E no , finding a 
fundamental system of solutions from the system: 

F* E a = 0, a e R fc , 

being: E c = { ve E n : v = E a, Va £ i? fc } 
and E no = { v G E" : F*v = } 

Being the fundamental system of solutions lineally idependent: {a\ a 2 ■ ■ ■ a q }, 
where: q < min ( k, n p ) , 
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dim E no = n - p = 1 

g i = J2 a ij e j, i=1 r • • ,9, E c n E no =( gt, ■ ■ ■ ,g k ) = E cno 

4- To complete the vectorial g i , i = 1,- • ■ ,q until a base in E c applying the 
method of elimination of Gauss of the matrice [ gi g2 ■ • • g q eie2 • • • Bk ] 

Note: When applying the elimination method all the vectors 
g i , i = 1,- • • , q. were used to build the escalar matrice. The vectors e^ 
j = 1,- • ■ ,k -q that are used in this submatriz is those that complete the 
base. 

5- To find a base in E no =(h\ h 2 ■ ■ ■ hi) 

6- To complete the vectorial g i , i = 1, • ■ • ,q until a base in E no, applying 
the one 

method of elimination of Gauss of the matrice [ gi g 2 • • • g q hih 2 ...h;] 
Note: All the vectorial g i , i = 1,- ■ • ,q was used, that form the main step, 

the vectors, h^. , j = 1,- • • , k-q that are used to complete this matrice is those 

that complete the base. 

7- To complete until a base of E n if it is necessary ( k + 1 - q < n ) ap- 
plying the elimination method of Gauss to the matrice [ gi g2 ■ • • g 9 ej 2 • • ■ 
e i k - q h n h n ■ ■ ■ hj l _ q v 1 v 2 ---v n ] 

where the v i , i = 1,- • • ,n are canonical vectors. To form the submatriz step 
only the v, were used. When it is not possible to use the remaining columns 
these Vj used are those that complete the base. 

8- Construct the matrice of change of base: 

If k + 1 - q = n then T =[ gi g 2 • • • g q e lx e l2 ■ ■ ■ e lk _ q h n h h ■ ■ ■ h n _ q ] 

Ifk + l-q< nthen T=[ gi g 2 • • • g q e ll e l2 ■■■ e lk _ q h jl h h ■■■h jl _ q v mi v m2 ■ ■ ■ 

9- To make the reduction to the canonical form finding 
^4=T -1 AT, B = T - x B y C = CT 

The state variables with the index { ii i 2 • • • ik-q } are completely control- 
lable and totally observables and they constitute the states of minimal repre- 
sentation of the analyzed system. 

In the canonical representation, we divide the state variables in 4 groups: 

Controllable and non observables, controllable and observables, not control- 
lable and non observables, not controllable and observables. 

To determine the efficiency of this algorithm you can analyze their temporary 
complexity. The number of operations is of order K 3 , where K = max (n r,n m), 
since in all the procedures you can apply the method of elimination of Gauss 
and the number of operations in this is of the order of the third power of the 
maximum between the number of rows and the number of columns [6] , for what 
its temporary complexity is of the order of K^, it is therefore an algorithm 
polynomial. 

Notice you the only jump in the logical structure of this exhibition takes 
place when they are accepted without demonstration, the practical approaches 
of control and observation total, really it is enough to demonstrate only one of 
these criterion, for the relationship of duality that exists among them. These 
criterions are of great practical utility for their algorithm form however they 
have an abstract character and the demonstrations are sophisticated points. 
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The demonstration of the necessity of the condition on the range of the 
control matrice can be demonstrated in a relatively simple way, if we lean on 
in the Theorem of Hamilton-Cayle [7] that in spite of being a very important 
theorem of the algebra, and specially useful inside the Theory of Control, and it 
is not included in the undergraduate courses. Several variants of demonstration 
of this Theorem appear in [8] , among them a very simple and elegant one based 
on purely algebraic arguments. 

The demonstration of the sufficiency of the criterion is more complex, a 
relatively simple variant can be found in [7] based on the concept of group of 
accessibility and using the theory of the differential equations. 

In [9] the necessity is demonstrated using Hamilton Cayley's Theorem and 
the formula of Sylvester that it is a generalization of the formula of interpolation 
of Lagrange for the case of the matrices. These information should be offered 
to the students, so that the interested ones in deepening in this theory have the 
appropriate bibliography. 

CONCLUSIONS 

It is developed an exhibition of the position and solution of the realization 
problem or identification of the lineal and stationary dynamic systems, without 
making use of the results of complex mathematical theories, related with the 
theory of the moments, the representation theorems of functional and lineal 
operators in spaces topologies and the theory of the matrices with elements in 
the field of the forms models - univaried that serve as foundation mathematics 
for the solution of this problem. 

The unit of the exhibition of the realization theory or identification is achieved 
through the intuition, based on the interpretation of such physical concepts as: 
overlapping, control and observation, the exemplification, the analysis, the syn- 
thesis and the comparison. 

A solution of the realization problem for the systems SISO is obtained in 
immediate form starting from the analysis and the comparison of the position 
and of the solution of the direct problem: " To find the transfer function 
applying the definition from this to the state pattern." 

For study of the systems MIMO study, one proceeds to decompose these, 
by virtue of the lineality, in a subsystem group SISO. each one of which are 
applied the previous results, obtaining a tena (A,B,C) in which thematrice A 
is a cellular diagonal matrice. The generalization of this result is expressed by 
means of the elaboration of an algorithm that allows to obtain (A,B,C) for any 
particular case. 

The existence of an infinite number of solutions of the equation matricial 
C (S I - A )- 1 B = G(S) 

puts on evidence in the fact that any system algebraically equivalent to the 
one obtained, also carries out the behavior entrance-exit defined by the transfer 
womb G(s). 

Starting from the theory of the subspaces of a vectorial space, the systems 
algebraically equivalent and the physical concepts of total observation and total 
control are based in the obtaining of the canonical decomposition of the obtained 



13 



pattern and is a completely controllable subsystem and totally observable and 
carries out the behavior pattern entrance-exit of G(s). 

The fact that the subsystem of the state variables totally observables and 
completely controllable it determines the realization minimal of G(s), it is easily 
accepted by the students, for the evident reduction of the dimension of the 
state space and the physical interpretation of the control concepts and total 
observation. 
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